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锡 林 河 流域 径流 量变 化 对 气候 变化 与 
人 类 活动 的 响应 
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摘 要 : 选取 锡 林 河 集 水 区 域 为 研究 区 ,基于 1968 一 2015 年 水 文 和 气象 数据 序列 , 采 
用 Mann-Kendall 突变 检验 和 双 累 积 曲 线 法 进行 突变 分 析 , 继 而 利用 双 累 积 曲 线 法 、 累 积 量 斜 率 变 化 


行 趋势 分 析 , 采 


治 区 河流 与 湖泊 生态 重点 实验 室 ,内 蒙古 
2. 蒙古 高 原生 态 学 与 资源 利用 教育 部 重点 实验 室 ,内 蒙古 
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率 比较 法 及 弹性 分 析 法 定量 研究 气候 变化 和 人 类 活动 对 径流 量 影响 的 贡献 率 。 人 研究 表明 :1968 一 2015 年 , 锡 林 河 


流域 年 径流 量 呈 显著 下 降 趋势 , 且 在 1984 年 和 
2000 年 和 2001 一 2015 年 则 作为 突变 期 ; 双 累 积 出 


2000 年 发 生 突变 。 基 于 此 ,将 1968—1984 年 
线 法 得 出 的 人 类 活动 对 径流 影响 贡献 率 介 于 86% ~88% ,与 其 他 
方法 所 得 结果 相差 较 大 ,而 累积 量变 化 率 比较 法 则 同时 考虑 降雨 和 蒸发 对 径流 的 影响 ,与 弹性 系数 法 得 出 的 结果 


FE 为 基准 期 1985 一 


基本 一 致 ,贡献 率 分 别 介 于 52.44% ~ 69. 02% (1985—2000 年 ) 和 42.39% ~43.64% (2001—2015 年 ) 。 综 合 而 言 ， 


不 同 定 量 方法 具有 不 同 的 基础 和 结构 ,使 
适 于 锡 林 河 流域 径流 量 的 定量 响应 研究 。 


包含 多 方法 集成 的 弹性 系数 法 相对 于 经 验 性 统计 法 而 言 更 可 靠 , 且 更 
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径流 量 对 气候 变化 与 人 类 活动 的 响应 是 现 阶段 
水 资源 演变 规律 研究 中 的 热点 科学 问题 5 。 目 前 ， 
开展 气候 变化 和 人 类 活动 对 径流 定量 响应 研究 方法 
主要 包括 经 验 统 计 法 .弹性 系数 法 和 水 文 模型 等 3 
类 中 。Ahn 等 中 选取 美国 4 个 区 域 为 研究 对 象 , 采 
用 线性 回归 分 析 法 水 量 平衡 法 水 文 模型 法 和 水 文 
敏感 性 分 析 法 4 种 方法 分 析 气 候 变 化 和 人 类 活动 对 
径流 的 影响 ,分 析 表 明 4 种 方法 在 4 个 研究 区 的 结 
果 接 近 ; Zeng 等 以 漳 河 流域 为 研究 区 ,从 年 、 季 、 
月 尺度 ,采用 模型 和 弹性 系数 法 研究 径流 量 对 气候 
变化 和 人 类 活动 的 定量 响应 ,结果 表明 湿 季 和 旱季 
季风 性 径流 量 对 人 类 活动 响应 的 贡献 率 分 别 为 


年 人 类 活动 对 于 径流 量 减 少 的 贡献 率 分 别 为 
73.70% 和 42.64% ; 焦 玮 研究 了 1963—2012 年 
锡 林 河 流域 河川 基 流 对 气候 变化 与 人 类 活动 的 响应 
特征 ,选取 1979 1991 年 和 1998 年 为 突变 年 份 , 划 
分 出 基准 期 和 措施 期 ,得 出 在 措施 期 内 人 类 活动 对 
基 流 的 贡献 率 分 别 为 45.85% .69.76% 和 56. 24% ; 
段 利 民 等 "研究 得 出 , 近 45 a 来 锡 林 河 流域 年 径 
流量 显著 下 降 ,1979 1986 年 和 1998 年 为 突变 年 
份 ,1980 一 1986 年 人 类 活动 对 径流 减少 的 贡献 率 为 
46.57% ;1987 一 1998 年 人 类 活动 对 径流 增加 的 贡 
献 率 高 达 91. 46% ;1999 一 2014 年 人 类 活动 对 径流 
减少 的 贡献 率 为 66.16% 。 尽 管 许多 学 者 已 在 锡 林 


57% 和 28% ;Wang 等 中 得 出 ,人 类 活动 对 黄河 径流 
变化 的 贡献 率 达 到 92. 07% ; 杨 立 哲 等 中 采用 M-K 
趋势 分 析 及 降雨 径流 双 累 积 曲 线 法 ,研究 了 锡 林 河 
流域 径流 量 的 年 内 、 年 际 变化 特点 ; 宋 小 园 等 ”2 针 
对 人 类 活动 和 气候 变化 影响 下 锡 林 河 的 水 文 响应 机 
制 , 选 择 基 准 期 为 1963—1985 年 ,在 不 考虑 潜在 组 
散 量 变化 的 影响 下 ,1986 一 195% 年 和 2000—2011 
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化 及 人 类 活动 对 流域 径流 量变 化 贡献 率 的 准确 估 
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图 1 锡 林 河 流域 示意 医 
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Fig.1 Schematic map of the Xilin River Basin 
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发 ""。 属 中 国 干 旱 、 半 干旱 地 区 典型 的 草原 型 河 
流 。 近 年 来 ,由 于 矿产 资源 过 度 开发 ,草原 生态 系统 
下 扑面 遭 到 严重 破坏 ,水 资源 时 空 分 布 格局 明显 改 
变 ,加 之 依赖 小 流量 河流 补给 的 草原 区 水 资源 系统 
的 脆弱 性 "” ,在 气候 变 暧 和 极端 事件 频 发 的 背景 
下 ,人 研究 锡 林 河 径流 变化 规律 ,定量 分 析 人 类 活动 和 
气候 变化 对 径流 量 的 贡献 率 , 有 助 于 了 解 内 陆 河流 
域 水 文 循环 过 程 的 演变 规律 。 本 文 在 已 有 人 研究 基础 
上 ,采用 经 验 统计 法 及 弹性 系数 法 两 大 类 7 种 方法 
开展 径流 量 对 气候 变化 和 人 类 活动 定量 响应 研究 ， 
以 期 使 研究 结果 更 为 准确 可 靠 。 


1 研究 区 概况 及 数据 来 源 


1.1 研究 区 概况 

锡 林 河 发 源 于 赤峰 市 克 什 克 腾 旗 宝 尔 图 西南 山 
顶 , 全 长 268. 1 km。 锡 林 河 流域 (图 1) 地 处 
115°32’ ~117° 16' E ,43°26' ~ 44"39'N ,流域 面积 
10 542 km? ,水 库 控制 面积 3 852 km? ,锡林浩特 水 
库 以 下 近 124.7 km 的 河流 已 断 流 。 根 据 1968 一 
2015 年 水 文 气 象 资料 统计 ,流域 以 锡林浩特 水 文 站 
日 径流 序列 实测 数据 作为 数据 源 ,多 年 平均 降水 量 


为 278.9 mm ,多 年 平均 蒸发 量 为 1 862.9 mm(20 
cm ZE) ;多 年 平均 气温 为 2.8 Y ,极端 最 高 气温 
为 31.6 C ,极端 最 低 气温 为 -33.1 °C ;多 年 平均 风 
速 为 3.4 m-s ,最 大 风速 为 14.5 m s; LFE 
均 径 流量 为 2.1 x10 ” m 。 伴 随 全 球 气候 变化 ,加 之 
过 度 放 牧 、 煤 矿 开 采 、 水 利 工程 建设 等 人 类 活动 , 锡 
林 河 流域 水 文 水 资源 问题 不 断 凸 显 ,河流 径流 量 显 
著 下 降 , 其 至 断 流 ,草场 退化 及 土壤 沙化 等 问题 亦 相 
继 频 现 。 
1.2 数据 来 源 

本 研究 水 文 数据 选取 锡林浩特 水 文 站 1968 一 
2015 年 的 日 径流 实测 数据 序列 ,气象 数据 选取 锡 林 
河流 域 气象 站 1968 一 2015 年 的 逐日 气象 数据 (包括 
降水 量 、 气 温 、 落 发 量 、 风 速 、 相 对 湿度 .气压 、 水 汽 
FR. A BRINE) 


2 研究 方法 


2.1 趋势 及 突变 分 析 

本 研究 运用 滑动 平均 法 、 累 积 距 平 法 ,线性 回归 
法 .Mann-Kendall 非 参数 统计 检验 法 ,揭示 水 文 气象 
数据 序列 的 趋势 性 变化 。 采 用 Mann-Kendall 突变 
检验 用 于 径流 突变 分 析 , 同时 应 用 双 累 积 曲 线 作 为 
突变 点 的 辅助 检验 。Mann-Kendall 检验 法 是 由 世界 
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张 艳 起 等 : 锡 林 河 流域 径流 量变 化 对 气候 变化 与 人 类 活动 的 响应 


气象 组 织 推荐 , Mann 和 Kendall 最 早 提出 -中 ,已 
成 功 应 用 于 气温 水质、 降水 和 径流 等 要 素 的 趋势 及 
突变 变化 分 析 中 "” 。 该 方法 的 特点 在 于 样本 无 需 
遵从 一 定 的 分 布 ,避免 少数 异常 值 的 干扰 ,在 水 文 、 
气象 等 非 正 态 分 布 的 数据 中 适用 性 良好 ,计算 简 
°° 。 

2.2 径流 量 对 气候 变化 和 人 类 活动 的 定量 响应 分 析 

对 于 一 个 确定 的 流域 ,径流 量 改 变 是 气候 变化 
和 人 类 活动 共同 作用 的 结果 ,可 表示 为 : 

AQ =AQc + AQ, (1) 
式 中 :AQ 为 径流 变化 总 量 ;AQc 为 径流 对 气候 变化 响 
应 的 变化 量 ;AQ; 为 径流 对 人 类 活动 响应 的 变化 量 。 
同时 ,通过 计算 突变 点 前 后 两 个 时 期 实测 径流 
量 的 差 值 , 即 可 得 到 径流 变化 总 量 , 即 : 

0=0 -0 (2) 
式 中 :01 为 基准 期 实测 径流 量 的 平均 值 ;Q, 为 突变 
期 实测 径流 量 的 平均 值 。 量 化 径流 量 对 气候 变化 和 
人 类 活动 的 响应 需要 计算 出 AQ: FAQ, 中 的 一 个 

因子 , 另 一 个 因子 即 可 通过 式 (1) 计 算得 出 。 

本 文 在 考虑 潜在 蒸 散发 和 实际 蒸 散发 的 基础 上 ， 

采用 2 种 经 验 统计 法 及 5 种 弹性 分 析 法 定量 分 析 气 
候 变化 和 人 类 活动 对 径流 量 的 影响 ,由 在 找 出 适 于 干 
旱 半 干旱 地 区 草原 型 河流 贡献 率 的 计算 方法 。 
2.2.1 传统 经 验 型 统计 法 “ 双 累 积 曲 线 法 是 目前 
应 用 于 水 文 要 素 一 致 性 或 长 序列 趋势 及 突变 分 析 中 
最 简单 直观 广泛 的 方法 ,通过 判定 两 个 变量 累 
积 值 之 间 直 线 斜 率 发 生 突变 的 位 置 , 即 可 确定 二 者 
RBA FAB IN Te] 。 

累积 量 斜率 变化 率 比 较 法 (SCRCQ 法 ) U” WA 
理 为 :基准 期 和 突变 期 累积 降水 量 一 年 份 线性 关系 
HRR Kp,» Kp, (mm + a’) ,斜率 变化 率 为 
Sp; 累积 蒸发 量 一 年 份 线性 方程 的 斜率 分 别 为 Ks、 


Kp, (mm a”) ,斜率 变化 率 为 5 ;累积 径流 量 一 年 
份 线性 方程 的 斜率 分 别 为 Kr, Kp, (mm + a`), R} 
率 变 化 率 为 54, 则 累积 变化 量 斜 率 变 化 率 为 突变 期 
累积 变化 量 一 年 份 线性 关系 的 斜率 差 值 与 基准 期 内 
累积 变化 量 一 年 份 线性 关系 的 斜率 比值 。 
以 降水 和 蒸发 代表 气候 变化 影响 的 变量 ,降水 
对 径流 量变 化 贡献 率 ( Cj ) 为 : 
_ 100 x 8, 


py (3) 
FRONT FE Tit Ba (EY BTR (Cy) A: 
100 x S, 
ET Sp (4) 
AŽ M SY AE EAER (Cy) A: 
C, =100 - C, - C; (5) 
上 述 式 中 :下 标 R 代表 径流 量 ;P 代表 降水 量 ;五 代 
表 蒸 发 量 。 
2.2.2 ”弹性 分 析 法 ”气候 变化 将 导致 降水 及 潜在 


蒜 发 能 力 的 变化 ,从 而 引起 径流 量 发 生 相 应 变化 。 
通过 分 析 径 流 对 降水 及 潜在 蒸 散 发 的 敏感 性 系数 ， 
径流 量 对 气候 变化 的 定量 响应 可 表示 为 : 


AQ, =8, SAP + ep, AE, (6) 
0 


SUH Ey 代表 潜在 蒸发 量 ;ep 和 er 分 别 为 径流 对 
于 降水 和 潜在 蒸发 的 弹性 系数 ;AP 和 AE 分 别 为 
降水 和 洪 在 蒸 散 发 的 变化 。 根 据 流域 长 时 段 水 量 平 
MARC =P -E,) 及 Budyko 假设 ,流域 实际 
蒸 散发 量 CE, ) 为 干 侣 指数 (由 = EAP) 的 函数 ,弹性 
系数 计算 公式 可 表示 为 : 


Hepte,,=1 (7) 


JEP Fb) I E (p) BYE LA F JLE Budyko 假设 
WH (#1). 


表 1 5 种 常用 的 基本 Budyko Rizh F(p) F'(p) RIAN 


Tab.1 Five basic Budyko assumptions in common use 


来 源 文献 F(o) Fg) 
21 gptanh (1/¢4) tanh(1/b) -4[ $b(e -ls +e”¢)?] 
_ F (tanh(1/¢) - seck? (1/6) /o 
22 an =e t S ; an. -e? 0:3 
[ dtan(1/p) (1 )] 0.5[ dtanh(1/p) (1 )] x ee 
1.5 
23 d+?) HG [1 +a% | } 
24 1+6-(1+¢°)'*,a=2.5 1- (14 gt) ¥@-' gen! 
25 (1 +wp)/(1 +wp+1/ġ),w=1 (w +2w/p-1 +1/ġ°)/(1 +w +1/)? 
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图 2 1968—2015 年 各 水 文 气象 要 素 年 际 趋势 变化 


Fig.2 Inter-annual trends of meteorological and hydrological factors during the period of 1968 - 2015 


3 结果 与 讨论 


3.1 水 文 气象 要 素 变化 特征 
3.1.1 水 文 气象 要 素 趋 势 分 析 ”从 图 2 和 表 2 可 
以 看 出 ,年 平均 气温 呈 增 加 趋势 ;年 平均 潜在 蒸发 量 
呈 显 著 上 升 趋势 ,线性 变化 速率 为 3.3 mm + a  , 变 
差 系数 为 0.05 , 变 幅 较 小 ;年 降水 量 和 汛期 降水 量 
变化 均 不 显著 ;年 径流 量变 化 呈 显 著 下 降 趋 势 ,21 
世纪 以 后 变化 极为 明显 。 
从 图 3 可 知 , 锡 林 河 流域 径流 量 呈 现 出 不 规则 
的 周期 性 波动 ,20 世纪 80 年 代 以 来 ,径流 量 波动 变 
化 比较 大 ; 20 世纪 80 年 代 后 期 和 90 年 代 前 期 径流 


R2 1968—2015 年 气象 和 水 文系 列 年 际 变化 及 
Mann-Kendall 趋势 检验 
Tab.2 Inter-annual variation of the meteorological and 
hydrological series during the period of 1968 -2015 and 
Mann-Kendall trend test 


气温 / 径流 量 降水 量 BERAE 
(C .a-!) /(mm + a7!) 
变化 趋势 0.04 -0.05 -0.03 3.30 
Zc 4.61 -3.39 -0.31 3. 80 
显著 水 平 ok ok kk S 六 六 六 


TE: * ** 表 示 在 0.01 水 平 上 显著 。 
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图 3 1968 一 2015 年 锡 林 河 流域 年 径流 与 降水 变化 过 程 
Fig.3 Variations of annual runoff volume and precipitation 


in the Xilin River Basin during the period of 1968 - 2015 


量 呈 现 增加 趋势 ;90 年 代 径 流量 高 于 多 年 平均 值 ; 
20 HR 21 世纪 初 径流 量 呈 现 极 显著 减少 趋势 。 
1968—2015 年 的 多 年 平均 径流 深 为 4. 19 mm ,最 大 
年 径流 深 10.25 mm(1993 年 ) ,最 小 年 径流 深 1. 68 
mm(2006 年 ) , 极 值 比 为 6.11。 

3.1.2 水 文 气象 要 素 突变 分 析 “应 用 M-K 突变 检 
验 法 分 析 1968 一 2015 年 锡 林 河 流域 年 气温 、 降 水 
量 潜在 蒸发 量 .径流 量 等 水 文 气 象 要 素 年 际 突 变 过 
程 , 给 定 显著 性 水 平 w =0.05, BI Up os = + 上 1.96 ,并 
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图 4 1968 一 2015 年 不 同 水 文 气象 要 素 年 际 M-K 突变 检验 


Fig.4 M-K mutation test of different meteorological and hydrological factors during the period of 1968 -2015 


利用 Anaconda 2.7 软件 得 到 UF UB ,结果 如 图 4 所 
示 。 由 图 4 可 以 看 出 ,1986 一 2015 年 锡 林 河 流域 气 
温 呈 显著 上 升 趋势 ,在 1988 年 发 生 突变 (图 4a) ; 降 
水 波动 不 显著 ,在 1979 年 出 现 突变 (网 4b) ;潜在 蒸 
发 呈 显 著 上 升 趋 势 ,在 1988 年 发 生 突变 (图 4c) ; 径 
流量 呈 显 车 减少 趋势 ,分 别 在 1984 年 和 2000 年 发 
生 突 变 (图 4d) ,这 与 宋 小 园 BEE? OSE E 
接近 ,20 世纪 90 年 代 出 现 一 定 程 度 的 增加 ,从 2001 
年 开始 又 出 现 明显 的 减少 趋势 ,上 且 超 过 0. 05 显著 性 
水 平 (mu = -1.96 ) ,甚至 超过 0. 01 显著 性 水 平 
(Usor = -2.58 ) ,表明 锡 林 河 流域 径流 量 减少 趋势 
极其 显著 。 其 中 当 显著 性 水 平分 别 设 定 为 0. 01、 
0.05 0.1 时 ,对 应 的 统计 量 临 界 检验 值 分 别 为 
+2.58,+1.96 和 + 上 1.65。 

总 体 而 言 , 锡 林 河 流域 气温 潜在 蒸发 量 均 呈 增 
加 趋势 ,降水 量 、 径 流量 则 呈 减 少 趋势 。 同 时 ,采用 
径流 一 降水 双 累 积 曲 线 法 识别 径流 序列 发 生 突变 的 
临界 年 份 , 互 为 验证 ,结果 如 图 5 所 示 。 

由 图 5 可 以 看 出 ,在 人 类 活动 和 气候 变化 的 双 
重 作用 下 , 双 累 积 曲 线 斜率 在 1984 .2000 年 发 生 明 
显 变化 , 且 3 个 阶段 线性 拟 合 趋势 线 斜 率 呈 现 先 增 
加 后 减少 的 趋势 , R 分 别 为 0. 998 7、0. 998 4, 
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Fig.5 Double cumulative curves of annual precipitation and 
runoff volume in the Xilin River Basin during the period of 


1968 -2015 


过 程 分 为 3 个 阶段 :1968 一 1984 年 .1985 一 2000 年 
和 2001—2015 年 ,其 中 1985 年 之 前 ,径流 量变 化 主 
要 受气 候 变 化 尤其 是 降水 量变 化 的 影响 ,1968 一 
1984 年 为 近似 天 然 径流 阶段 ,看 作 基准 期 ,人 类 活 
动 对 锡 林 河 流域 的 干预 较 少 。 此 后 ,人 类 活动 的 影 
响 逐 渐 增 大 ,径流 量 的 变化 除了 受到 降水 量 等 自然 
因素 的 影响 外 ,还 包括 人 类 活动 的 影响 ;1985 一 2000 
年 .2001 一 2015 年 为 气候 变化 和 人 类 活动 双重 作用 
影响 变 强 阶 段 ,降水 对 径流 产生 影响 ,人 类 亦 通过 改 


0.994 2 ,尤其 是 第 一 阶段 (1968 一 1984 年 ) 和 第 二 
阶段 (1985 一 2000 4E) R? 均 超 过 0. 998 ,说 明 线 性 拟 
合 程度 较 高 。 综 上 所 述 ,本 研究 将 降水 一 径流 变化 


变 下 垫 面 间接 产生 影响 ,尤其 是 21 地 纪 以 来 径流 变 
化 显著 ,该 结果 与 M-K 检验 UF-UB 统计 结果 一 致 。 
由 图 6 可 见 , 随 着 时 间 推 移 , 1985 一 2000 年 年 
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Fig.6 Scatter diagram of annual runoff volume and precipitation 


in different periods 


径流 量 和 降水 量 二 者 关系 整体 向 右 推 移 ,说 明 相同 
量 级 降水 量 下 径流 量 相应 增加 ;然而 2001—2015 年 
二 者 关系 点 向 左 推移 且 极 其 显著 ,说 明 相同 量 级 降 
水 量 下 径流 量 明显 减少 ,间接 验证 了 非 降水 因素 对 


径流 量 影响 的 可 能 性 。 


3.2 径流 量 对 气候 变化 和 人 类 活动 的 定量 响应 
3.2.1 双 累 积 曲线 法 定量 区 分 径流 变化 影响 因子 


根据 M-K 突变 检验 结果 ,将 研究 时 段 划分 为 3 个 


阶段 :1968 一 1984 年 (基准 期 ) , 1985—2000 年 ( 38 
变 期 ) .2001 一 2015 年 (突变 期 ) 。 

年 平均 降水 序列 与 年 平均 径流 量 的 线性 回归 方 
程 可 表示 为 : 


R=0.014 6P -0.3444 (8) 


锡 林 河 流域 不 同时 有 段 理论 年 均 径流 量 可 以 由 式 
(8) 计 算 求 得 ,突变 期 实测 径流 量 与 基准 期 实测 值 
之 间 的 差 值 包括 人 类 活动 影响 和 降水 变化 影响 两 部 
分 ,其 中 突变 期 降水 量 对 径流 量变 化 的 影响 值 可 由 
基准 期 实测 值 与 各 时 段 的 计算 值 差 值得 到 。 计 算 结 
果 见 表 3。 
3.2.2 累积 量 斜 率 变 化 率 分 析 方法 区 分 径流 变化 
影响 因子 “该 方法 中 ,以 1984 .2000 年 2 个 突变 年 
份 为 界 , 将 各 变量 累积 量 的 变化 阶段 划分 为 3 个 阶 
BE, 1968—1984 年 为 Ta, 1985—2000 年 为 7b, 
2001—2015 年 为 Te ,根据 SCRCQ 法 计算 各 阶段 线 
性 关系 斜率 的 变化 率 ( 表 4) 。 根 据 公式 分 别 估算 各 
响 因 素 对 径流 量变 化 的 贡献 率 ,结果 见 表 5。 
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RI 锡 林 河 流域 降水 和 人 类 活动 对 径流 影响 的 计算 结果 


Tab.3 Calculated results of the effects of precipitation and human activities on runoff volume in the Xilin River Basin 


年 均 径 流量 /mm ; 降水 因子 人 类 活动 因子 
起 止 年 份 。 降水 量 /nm © te 一 ee 
实测 值 计算 值 影响 值 影响 率 /% 影响 值 影响 率 /% 
1968—1984 282.08 4.45 
1985—2000 289. 36 5.47 4.57 -1.02 -0.12 11.67 -0.90 88. 33 
2001—2015 264.12 2.54 4.20 1.91 0.25 13.09 1.66 86.91 
表 4 不 同时 段 各 累积 量 斜率 变化 率 
Tab.4 Change rates of cumulative slope in different periods 
时 期 径流 量 R 降水 量 P RAEE E 
B RFE Kp 变化 率 Sk/% 斜率 Kp 变化 率 Sp/% 斜率 Kg 变化 率 Sg/% 
Ta 4.64 296. 18 1 006.9 
Tb 5.79 24. 82 302. 49 2. 13 1 104. 3 9. 67 
Tc 2.46 -57.56 259. 89 -14.08 901.58 -18.36 
表 5 不 同时 段 气 候 变化 和 人 类 活动 对 径流 影响 的 贡献 率 。” 子 对 径流 变化 的 影响 ,1985 一 2000 年 的 计算 结果 见 


Tab.5 Contribution rates of climate change and human 


activities to runoff volume in different periods 


时 期 Cp + Cg/% Cp/% 

Th 47. 56 52. 44 

Te 56. 36 43.64 
3.2.3 弹性 系数 法 定量 区 分 径流 变化 影响 因子 


按照 年 时 间 尺 度 , 采 用 5 种 基于 Budyko 假设 的 函数 
形式 ,计算 不 同 阶段 径流 对 降水 的 弹性 系数 及 各 因 


表 6,2001 一 2015 年 的 计算 结果 见 表 7。 

3.3 讨论 

3.3.1 不 同方 法 比较 ” 表 8 对 比 了 双 累 积 曲 线 、 
累积 量变 化 率 比 较 法 及 基于 弹性 系数 法 (5 种 方法 
平均 值 ) 的 计算 结果 。 对 双 累 积 曲 线 、 累 积 量 斜率 
变化 率 分 析 和 弹性 系数 3 种 分 析 方 法 ,分割 出 人 类 
活动 影响 径流 量 均 采 用 年 径流 量 和 年 降水 总 量 。 在 
2 个 突变 期 内 , 双 累 积 曲 线 法 得 出 的 人 类 活动 对 径 
流 影响 贡献 率 介 于 86% ~88% , 与 其 他 方法 所 得 结 
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表 6 1985 一 2000 年 弹性 系数 法 对 气候 变化 和 人 类 活动 
影响 的 定量 区 分 

Tab.6 Quantitative distinction between climate change 

and human activities during the period of 1985 —2000, 


based on the elastic coefficient method 


ep AQo/mm AQy/mm AQe/% AQn/% 
Oldekop 2.98 0.38 0.64 37.06 62.94 
Budyko 3.75 0.47 0.54 46.63 53.37 
Pike 2.96 0.37 0.64 36.85 63.15 
(HBR 2.49 0.32 0.70 30.98 69.02 
Zhang 等 2.82 (0.36 0. 66 35.12 64.88 
平均 值 ”3.00 ”0.38 0. 64 37.33 62.67 


表 7 2001—2015 年 弹性 系数 法 对 气候 变化 和 人 类 活动 
影响 的 定量 区 分 
Tab.7 Quantitative distinction of the effects of climate 
change and human activities during the period of 


2001 -2015 based on the elastic coefficient method 


ep AQc/mm AQp/mm AQce/W = AQy/% 


Ol’dekop 2.98 -1.10 -0.81 57.61 42.39 


Budyko 3.63 -1.38 -0.53 72.43 27.57 

Pike 2.97 -1.09 -0.82 57.27 42.73 
(ELBE 2.49 一 0.89 一 1.02 46.41 53.59 
Zhang 等 2.84 一 1.04 -0.87 54.27 45.73 
平均 值 2.98 -1.10 -0.81 57.60 42.40 


表 8 不 同 阶段 3 种 方法 对 气候 变化 和 人 类 活动 
影响 的 定量 区 分 
Tab.8 Quantitative distinction of the effects of climate 
change and human activities at different stages based on 
three methods 
累积 量变 化 率 
比较 分 析 法 


气候 变 ATR ”气候 变 人 类 活 ARE 人 类 活 
化 /% 动 /% r% 动 /% 化 /% 动 /% 


1985—2000 11.67 88.33 37.33 62.67 47.56 52.44 


双 累 积 曲线 弹性 系数 法 


时 Be 


2001—2015 13.09 86.91 57.60 42.40 56.36 43.64 


果 相 差 较 大 ,原因 在 于 该 法 仅 考虑 了 降水 对 径流 的 
影响 ,并 未 考虑 蒸发 量 等 其 他 因素 的 影响 ; 而 累积 量 
变化 率 比 较 法 则 同时 考虑 降雨 和 蒸发 对 径流 的 影 
响 , 如 果 量 化 出 该 法 中 降水 对 径流 的 贡献 率 
(8.58% All 24.46% ) ,其 与 双 累 积 曲 线 中 降水 对 径 
流 贡献 率 (11.67% 和 13.06% ) 相 近 。 此 外 ,弹性 系 
数 法 和 累积 量变 化 率 比较 法 得 出 的 结果 基本 一 致 ， 


这 些 方法 计算 得 到 的 人 类 活动 对 径流 影响 的 贡献 率 
分 别 为 52. 44% ~ 69. 02% (1985 一 2000 年 ) 和 
42.39% ~43. 64% (2001 一 2015 年 ) 。 与 经 验 性 统 
计 法 相 比 ,弹性 系数 法 不 仅 考 虑 了 潜在 蒸发 的 影响 ， 
且 基 于 Budyko 水 热平衡 耦 合理 论 和 水 量 平 衡 方程 
分 析 径 流 对 于 降水 和 蒸发 的 敏感 性 ,与 真实 水 循环 
过 程 更 接近 。 然 而 ,基于 不 同 的 Budyko 假设 函数 也 
可 能 得 出 不 甚 一 致 的 结论 ,如 Ol’ dekop , (Hu FX , 
Zhang “ .Ture-Pike 与 Budyko 得 出 的 结论 之 间 的 差 
别 ,原因 在 于 基于 不 同 Budyko 水 热平衡 耦合 理论 所 
包含 的 函数 参数 变化 。 基 于 此 ,采用 平均 值 可 减少 
不 确定 性 ,用 于 定量 评估 更 加 合理 。Wu 等 “采用 
多 种 方法 定量 分 析 气 候 变化 和 人 类 活动 对 径流 的 影 
响 ,同样 得 出 弹性 系数 法 所 得 到 的 结果 更 加 合理 。 
综 上 所 述 , 不 同 定量 方法 具有 不 同 的 基础 和 结构 , 采 
用 包含 多 方法 集成 的 弹性 系数 法 相对 于 经 验 性 统计 
法 而 言 ,更 加 准确 可 靠 。 

3.3.2 不 同 阶段 人 类 活动 及 气候 变化 对 径流 量 贡 
献 率 变化 原因 分 析 “由 以 上 分 析 可 知 ,1984 一 2000 
年 锡 林 河 流域 人 类 活动 影响 大 于 气候 变化 ,2000 年 
以 后 人 类 活动 影响 呈 小 幅 降 低 趋 势 ,2001 一 2015 年 
气候 变化 对 径流 影响 略 大 于 人 类 活动 ,为 气候 变化 
和 人 类 活动 双重 作用 变 强 的 阶段 ,该 结果 与 焦 
HO EO 等 的 研究 结果 基本 一 致 。 

1985—2000 年 径流 量 呈 现 一 定 程 度 的 增加 趋 
势 ,与 基准 期 相 比 增加 22. 84% 。 其 原因 可 能 是 : 自 
20 世纪 80 年 代 中 期 开始 ,人 口 增 长 速度 有 所 减缓 ; 
相 较 1968—1984 年 及 2001—2015 年 2 个 时 段 ,该 
时 段 锡林郭勒 盟 牲 畜 总 量 增加 了 将 近 2 倍 ,过 度 放 
牧 .超载 养 冀 ,水 利 工程 等 人 类 活动 导致 锡 林 河 流域 
下 热 面 改 变 , 削 弱 了 下 执 面 对 雨水 的 拦蓄 作用 ， 
在 一 定 程度 上 会 影响 径流 规律 ,造成 径流 量 加 大 。 
1985 一 2000 年 流域 降水 增加 ,造成 锡 林 河 流域 径流 
量 增 加 3 600 m ,潜在 蒸发 量 的 减少 造成 锡 林 河 流 
域 径 流量 增加 3 550 m ,人 类 活动 造成 流域 径流 量 
增加 10 m°, 

2001 一 2015 年 径流 量 呈 现 明 显 的 减少 趋势 , FE 
变 点 明显 ,与 基准 期 相 比 降幅 达到 42.4% ,气候 变 
化 影响 略 大 于 人 类 活动 ,说 明 人 类 活动 和 气候 变化 
的 双重 作用 导致 径流 量 的 减少 。 就 气候 变化 而 言 ， 
气温 突变 和 趋势 检验 结果 表明 ,21 世纪 以 来 气温 呈 
明显 上 升 趋势 , 与 全 球 变 暖 趋势 相 一 致富 ,加 之 降 
水 量 减少 6.37% ,潜在 蒸发 量 增加 5. 88% ,导致 径 


流量 在 一 定 程度 上 减少 。 就 人 类 活动 而 言 ,尽管 从 
2000 年 开始 ,锡林郭勒 盟 开展 了 一 系列 的 生态 保护 
工作 ,包括 “ 双 权 一 制 "“ “退耕 还 林 还 草 ”“ 围 封禁 
牧 ”“ 天 然 林 保 护 工程 "“ 京 津 风沙 源 治理 ”等 水 
土 保持 措施 ,防治 了 水 土 的 流失 ,增加 了 鞭 水 ,一定 
程度 上 减轻 了 人 类 活动 对 径流 的 影响 ;还 通过 增加 
植被 盖 度 或 林 草 生长 ,减少 出 口 断面 的 径流 量 。 然 
而 ,第 二 三 产业 的 迅 独 发 展 , 尤 其 是 煤 训 .石油 等 矿 
产 资源 的 过 度 开采 , 矿 点 较 多 从 而 消耗 水 量 , 地 表 水 
利用 量 及 地 下 水 开采 量 不 断 加 大 ,加 之 灌溉 引水 强 
度 的 增加 ,导致 2000 年 后 锡 林 河 流域 径流 量 明显 下 
降 。 根 据 基 于 Budyko 水 热平衡 耦合 理论 的 弹性 系 
数 法 量化 分 析出 降水 及 莹 发 对 径流 的 影响 ,从 而 量 
化 出 人 类 活动 对 径流 量 的 影响 。2001 一 2015 年 流 
域 降 水 减少 和 潜在 蒸发 量 增 加 ,造成 锡 林 河 流域 径 
流量 减少 1. 04 x 107 m ,人 类 活动 造成 流域 径流 量 
减少 8 700 mm。 


4 ”结论 与 展望 


结论 
(1) 锡 林 河 流域 气温 ,潜在 蒸发 均 呈 显著 增加 
趋势 ,降水 量 波动 不 明显 ,径流 量 呈 显著 减少 趋势 。 

(2) 径流 量 在 1984 年 和 2000 年 发 生 突变 , 基 
于 此 ,将 1968—1984 年 作为 基准 期 ,1985 一 2000 年 
和 2001 一 2015 年 作为 突变 期 。 其 中 , 1985 一 2000 
年 径流 呈 增 加 趋势 , 较 基 准 期 增加 22. 84% ,同时 降 
水 增加 5.2% , 游 在 蒸发 量 减 少 1. 4% ;2001 一 2015 
年 径流 呈 减 少 趋势 , 较 基 准 期 减少 28.9% ,同时 降 
水 增加 2.58% ,潜在 蒸发 量 减少 0.01% 。 

(3) 突变 期 1985—2000 年 ,利用 双 累 积 曲 线 
法 .累积 量 斜率 变化 率 分 析 法 及 弹性 系数 法 计算 的 
人 类 活动 对 径流 量 影响 的 贡献 率 分 别 为 88.33% 、 
52.44% 和 62. 67% ,人 类 活动 在 该 时 期 径流 变化 中 
起 主导 作用 ,造成 流域 径流 量 增加 6 600 m ;突变 期 
2001 一 2015 年 ,3 种 方法 计算 的 人 类 活动 对 径流 量 
影响 的 贡献 率 分 别 为 86. 91% .43. 64% 和 42.40% , 
双 累 积 曲线 法 数值 偏 大 ,累积 量 斜 率 变 化 率 分 析 法 
与 弹性 系数 法 得 出 的 结论 基本 一 致 ,气候 变化 较 人 
类 活动 对 径流 的 影响 略 大 ,人 类 活动 约 造 成 流域 径 
流量 减少 8 700 m 。 依 据 方法 的 基础 和 结构 ,采用 
包含 多 方法 集成 并 配合 水 热平衡 耦合 理论 的 弹性 系 
数 法 相对 于 经 验 性 统计 法 而 言 ,更 加 准确 可 靠 。 


4.1 
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4.2 展望 

(1) 本 文采 用 经 验 统计 法 和 弹性 系数 法 开展 径 
流量 对 气候 变化 与 人 类 活动 的 定量 响应 分 析 , 其 中 
所 涉及 的 水 文 气象 参数 有 限 , 未 来 可 使 用 分 布 式 或 
半分 布 式 水 文 模型 进一步 开展 相应 研究 ,以 期 得 到 
更 为 准确 可 靠 的 研究 结论 。 

(2) 人 类 活动 影响 对 径流 量变 化 不 容 忽 视 , 但 
由 于 社会 经 济 资料 的 缺乏 ,文中 未 深入 定量 揭示 已 
实施 的 生态 保护 水土 保持 , 亦 或 资源 过 度 开发 等 对 
径流 量 造成 的 影响 ,并 探讨 其 机 制 , 未 来 可 进一步 深 
化 该 领域 的 研究 ,以 期 揭示 人 类 活动 对 径流 量 影响 
的 症结 所 在 ,为 流域 生态 修复 提供 决策 依据 。 
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Response of Runoff Volume Change to Climate Change and 
Human Activities in the Xilin River Basin 


ZHANG Yan-xia', YU Rui-hong’*?, XUE Hao’, DUO Lan', HU Hai-zhu'*, LV Xi-xi’ 
(1. College of Ecology and Environment/Inner Mongolia Key Laboratory of River and Lake Ecology , 
Inner Mongolia University , Hohhot 010021 ,/nner Mongolia , China; 
2. Key Laboratory of Mongolian Plateau Ecology and Resources Utilization under Ministry of Education , 
Hohhot 010021 Inner Mongolia , China ) 


Abstract: Study on the quantitative response of runoff volume to climate change and human activities can not only 
provide a theoretical basis for the redistribution of water resources in a drainage basin, but also a policy reference for 
the management of water resources. This study was based on the hydrological and meteorological data series in the 
Xilin River Basin ,a typical grassland basin, during the period from 1968 to 2015 ,and the Mann-Kendall trend test 
was used to analyze the change trend of runoff volume. The Mann-Kendall mutation test and the double-accumula- 
tion curve method were also used. The contribution rates of climate change and human activities to the change of 
runoff volume were quantified using the double-accumulation curve method ,cumulant slope rate comparison method 
and elasticity analysis method. The results showed that the annual runoff volume in the Xilin River Basin was in a 
significant decrease trend during the period from 1968 to 2015 , and its mutations occurred in 1984 and 2000 respec- 
tively. On which the period of 1968 - 1984 was regarded as the reference period, and the periods of 1985 - 2000 
and 2001 -2015 were the mutation periods. The results from the double-accumulation curve method revealed that 
the contribution rate of human activities to runoff volume varied in a range of 86% - 88% ,which was quite different 
from the results obtained by other methods. The effects of both rainfall and evaporation on runoff volume were simul- 
tanously considered in the cumulant change rate comparison ,so the results estimated by which and by the elasticity 
coefficient method were similar, and the contribution rates varied in ranges of 52.44% -69.02% (1985 - 2000) 
and 42.39% -43.64% (2001 -2015) respectively. Holistically ,the foundations and structures of different quanti- 
tative methods were different ,the multiple methods were integrated in the elastic coefficient method, and this method 
was more reliable in quantifying the response of runoff volume in the Xilin River Basin. 


Key words: runoff volume; human activity; climate change; quantitative response; Xilin River Basin 


